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Abstract:  Landscape  assessment  of  soil  moisture  is  critical  to  understanding  the 
hydrological  cycle  at  the  regional  scale  and  in  broad-scale  studies  of  biophysi¬ 
cal  processes  affected  by  global  climate  changes  in  temperature  and  precipitation. 
Traditional  efforts  to  measure  soil  moisture  have  been  principally  restricted  to  in 
situ  measurements,  so  remote  sensing  techniques  are  often  employed.  Hyperspectral 
sensors  with  finer  spatial  resolution  and  narrow  band  widths  may  offer  an  alternative 
to  traditional  multispectral  analysis  of  soil  moisture,  particularly  in  landscapes  w  ith 
high  spatial  heterogeneity  .  This  preliminary  research  evaluates  the  ability  of  remotely 
sensed  hyperspectral  data  to  quantity  soil  moisture  for  the  Little  River  Experimental 
Watershed  (LREW),  Georgia.  An  airborne  hyperspectral  instrument  with  a  short- 
wavelength  infrared  (SWIR)  sensor  was  flown  in  2005  and  2007  and  the  results  were 
correlated  to  in  situ  soil  moisture  values.  A  significant  statistical  correlation  ( R 2  value 
above  0.7  for  both  sampling  dates)  for  the  hyperspectral  instrument  data  and  the  soil 
moisture  probe  data  at  5.08  cm  (2  inches)  was  determined.  WTiile  models  for  the 
20.32  cm  (8  inches)  and  30.48  cm  (1 2  inches)  depths  were  tested,  they  were  not  able 
to  estimate  soil  moisture  to  the  same  degree. 


INTRODUCTION 

Soil  moisture  is  a  critical  process  in  the  water  cycle  and  its  assessment  is  of  para¬ 
mount  importance  in  forecasting  changes  in  the  water  balance  of  a  region  (Salvucci 
et  al.,  2002).  In  agricultural  production,  the  spatial  variability  of  soil  moisture  can  be 
responsible  for  low  or  spatially  variable  crop  yields,  as  soil  moisture  is  required  to 
make  nutrients  soluble  for  plant  absorption. 

Soil  moisture  fluctuates  both  spatially  and  temporally  due  to  factors  such  as  soil 
type,  soil  horizon,  and  other  site-specific  geologic  and  climatic  conditions.  Traditional 
efforts  to  measure  soil  moisture  have  been  principally  restricted  to  in  situ  measurements 
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(Engman,  1999)  and  predominately  reliant  on  direct  field  instrumentation  and  point 
data  (Steele-Dunne  et  al.,  2010),  which  can  be  costly  when  used  to  obtain  spatially 
distributed  data  over  a  large  region  or  landscape  (Jensen  and  Hodgson,  2004).  As  an 
alternative  to  ground  point  data,  remote  sensing  data  are  seen  as  a  promising  approach 
to  evaluate  soil  moisture  characteristics  in  different  landscapes  and  sampling  areas  by 
providing  in  a  regional  description  of  the  water  redistribution  at  different  temporal  and 
spatial  resolutions  (McCabe  and  Wood,  2006). 

The  advantage  of  remote  digital  imaging  is  that  relatively  frequent  measurements 
can  be  made.  By  subtracting  spatially  referenced  data,  time-displacement  associations 
(i.e.,  time  series)  can  be  produced.  Beginning  in  the  1920s,  aerial  photographs  were 
used  to  delineate  and  map  boundaries  of  different  soil  series  by  identifying  dissimilari¬ 
ties  in  photographs  (Helms  et  al.,  2002).  By  focusing  primarily  on  using  vegetation 
as  a  proxy  for  soil  water  holding  capacity,  mapped  soils,  derived  from  those  image 
demarcations,  can  aid  managers  and  farmers  in  identifying  homogenous  zones  within 
agricultural  fields  for  guiding  nutrient  application  rates  using  variable  rate  technology' 
(VRT)  (Baumgardner  et  al.,  1985;  Walthall  et  al.,  2001;  Gish  et  al.,  2002;  Daughtry'  et 
al.,  2002). 

More  recently,  soil  scientists  have  used  multispectral  images  for  mapping  land 
cover  attributes  and  related  moisture  conditions.  The  wetness  in  the  top  few  centime¬ 
ters  of  soil  can  be  measured  remotely  using  multispectral  remote  sensing  or  certain 
radar  bands  (Njoku  and  Entekhabi,  1 996).  Transformations  of  the  spectral  reflectance 
in  remotely  sensed  images  can  provide  useful  information  on  soil  water  content  and, 
if  augmented  with  existing  soil  and  other  geospatial  information,  can  provide  spatially 
representative  estimates  of  soil  moisture.  Most  studies  using  remotely  sensed  data  to 
look  for  relationships  with  soil  properties  have  focused  on  the  0.3  to  2.8  micrometer 
region  of  the  electromagnetic  spectrum  (Bushnell,  1932;  Barnes  et  al.,  2003). 

Despite  numerous  studies,  a  universal  spectral  range  and  sampling  scheme  has 
yet  to  be  determined  for  best  detecting  soil  moisture.  Visible  near-infrared  (VNIR) 
(0.4-1. 4  pm),  near-infrared  (NIR)  (0.75-1.4  pm),  short-wavelength  infrared  (SW1R) 
(1. 4-3.0  pm),  thermal  infrared  (3.5-20  pm),  and  microwave  (0.3-300  cm)  have  each 
shown  some  promise.  In  general,  the  literature  suggests  that  the  longer  the  wavelength 
the  better  the  potential  ability  to  predict  soil  moisture.  NIR  appears  to  be  a  suitable 
option  for  bare-soil  fields  after  a  rainfall  (Milfred  and  Kiefer,  1976)  if  the  particle  size 
and  soil  type  are  the  same  as  the  calibration  set  (Slaughter  et  al.,  2001).  NIR,  when 
used  in  conjunction  with  visible  bands,  holds  promise  to  categorize  soil  samples  by 
moisture  content,  especially  when  soil  variability  is  low  (Mouazen  et  al.,  2006).  As 
a  better  alternative,  SWIR  has  been  shown  to  be  more  effective  than  VNIR  at  moni¬ 
toring  soil  moisture  changes  (important  to  predicting  soil  moisture),  as  long  as  the 
soil  moisture  values  remain  at  less  than  or  equal  to  50%  of  volumetric  water  content 
(Lobell  and  Asner,  2002),  especially  the  creation  of  Gaussian  models  for  soils  from  the 
Mediterranean  region  (Whiting  et  al.,  2004).  Measurements  from  the  thermal  infrared 
portion  of  the  electromagnetic  spectrum  have  also  been  conducted.  Thermal  infrared 
emissivity  increases  as  soil  moisture  content  decreases  at  low  water  content  levels 
(Mira  et  al.,  20 1 0)  and  thermal  infrared  images  with  high  spatial  resolution  can  be  used 
to  map  soil  moisture  at  a  fine  scale  (Katra  et  al.,  2006). 

Microwave  bands  also  have  the  capacity  to  measure  soil  moisture  (Moran  et 
al.,  1998;  Barnes  et  al.,  2003)  over  a  selection  of  topographic  and  vegetation  cover 
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settings  (Jackson  et  al.,  1997).  On  relatively  unvegetated  soil  it  is  possible  to  obtain 
reasonably  accurate  soil  moisture  estimates  using  active  microwave  X-  (2.4-3. 8  cm) 
and  C-band  (3. 8-7.5  cm)  radar  imagery  (Jensen  et  al.,  2005).  Because  of  the  signal 
ability  to  penetrate  the  earth,  microwave  technology  allows  scientists  to  gather  data 
on  the  properties  of  the  zone  from  the  soil  surface  to  the  groundwater  table  (Crow  et 
al.,  2005).  Extending  this  technology  to  routine  measurements  from  a  satellite  system 
would  supply  notable  spatial  and  temporal  characteristics  to  hydrologic  measurements 
(Engman,  1999).  Additional  microwave  remote  sensing  investigations  specific  to  our 
study  area  are  expanded  upon  in  the  next  section  of  this  paper. 

As  an  alternative  to  broad-band  multispectral  remote  sensing  sensors,  which 
collect  spectral  information  across  relatively  wide  ranges  of  the  electromagnetic  spec¬ 
trum,  hyperspectral  sensors  collect  remote  sensing  data  as  images,  with  each  image 
representing  a  different  spectral  band  with  higher  spectral  resolution  than  broad-band 
sensors.  When  stacked  together,  the  two-dimensional  spatial  (image)  data  sets  create 
hyperspectral  data  cubes,  with  two  dimensions  containing  spatial  information  and 
the  third  dimension  containing  spectral  information  (Warren  et  al.,  1997;  Vogt  et  al., 
2004). 

The  possibilities  of  using  hyperspectral  data  to  obtain  soil  water  content  are  still 
uncertain,  but  appear  promising.  For  example,  Grandjean  et  al.,  (2010)  identified  air¬ 
borne  hyperspectral  sensing  as  one  of  the  anticipated  geophysical  sensor  technologies 
that  will  aid  the  dilemma  of  performing  soil  data  collections  at  the  catchment  scale.  The 
use  of  airborne  hyperspectral  sensing  will  bridge  a  technological  gap  for  the  purpose 
of  developing  steadfast  and  cost-effective  geophysical  mapping  solutions,  especially 
for  data  collection  in  digital  soil  property  mapping  (soil  texture,  soil  water  content,  soil 
hydraulic  properties,  bulk  density,  and  soil  organic  matter).  While  airplane-mounted 
sensors  have  distinct  advantages  over  satellite-based  systems,  the  approach  has  yet  to 
be  thoroughly  tested.  The  ability  to  have  data  with  finer  spatial  resolution  and  more 
control  over  temporal  resolution  from  aircraft  mounted  sensors  is  an  advantage  over 
satellite-based  data.  The  larger  ground  area  covered  in  less  time  is  a  benefit  over  in 
situ  measurements. 


OBJECTIVES 

The  purpose  of  this  research  was  to  test  the  statistical  correlation  between  spectral 
bands  of  an  airborne  hy  perspectral  sensor  and  soil  moisture  ground  readings  at  three 
different  depths.  The  research  explored  whether  the  reflectance  values  available  in 
hyperspectral  remote  sensing  datasets  can  be  used  in  the  monitoring  of  soil  water  con¬ 
tent.  This  study  contributes  to  an  ongoing  effort  to  obtain  spatially  accurate  soil  mois¬ 
ture  data  using  remote  sensors  that  can  be  used  in  hydrological  and  climatic  research 
on  such  issues  as  drought  prevention/risk,  analyzing  ground-atmosphere  feedback, 
and  crop  irrigation  techniques  among  others. 

STUDY  AREA 

The  Southeast  Watershed  Research  laboratory  (SEWRL)  of  the  U.S.  Department 
of  Agriculture  has  30  instrumented  stations  continuously  collecting  soil  water  con¬ 
tent  for  the  Little  River  Experimental  Watershed  (LREW),  a  small  watershed  in  south 
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Fig.  1.  L  .ittle  River  Experimental  Watershed  (LREW)  and  the  Suwannee  River  Basin  of  southern 
Georgia  and  northern  Florida. 


Georgia  (Bosch  et  ah,  2007).  The  LREW  is  a  sub-basin  of  the  Suwannee  River  Basin 
of  southern  Georgia  and  northern  Florida  (Fig.  1).  The  study  area  has  a  mild  climate 
with  warm,  humid  summers  and  mild  winters,  and  stream  flow  responses  that  are 
dominated  by  annual  and  seasonal  precipitation  trends  (NOAA,  2002;  Priest,  2004; 
Bosch,  2006b). 

Prior  studies  conducted  on  sample  plots  in  this  area  show  high  soil  moisture  vari¬ 
ability  between  plots  (Giraldo  et  al.,  2008).  However,  within  plots  there  is  high  homo¬ 
geneity  and  temporal  stability  when  using  soil  moisture  readings.  In  addition,  soil 
moisture  measurements  vary  with  different  land  use/land  cover  types  (Giraldo  et  al., 
2008,  2009). 

PREVIOUS  MICROWAVE  INVESTIGATIONS  CONCERNING  SOIL 
MOISTURE  IN  THE  LITTLE  RIVER  EXPERIMENTAL  WATERSHED 

(LREW) 

As  previously  stated,  Jensen  et  al.  (2005)  and  Crow  et  al.  (2005)  have  shown 
that  soil  moisture  estimates  using  microwave  remote  sensing  techniques  are  possible, 
especially  at  the  soil  surface,  but  they  can  also  be  expanded  to  a  few  centimeters  depth. 
Both  active  and  passive  sensors  supply  information  on  surface  reflectivity  that  in  turn 
is  related  to  soil  moisture  through  the  Fresnel  reflection  equations  (Ulaby  et  al.,  1986; 
Jackson,  2003).  Fresnel  reflection  equations  allow  microwave  sensor  measurements  to 
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be  equated  to  soil  moisture,  based  on  relationships  between  reflected  and  transmitted 
electromagnetic  waves. 

Data  were  collected  on  the  LREW  using  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  Polarimetric  Scanning  Radiometer  (PSR/CX)  sensor  during 
Soil  Moisture  Experiment  2003  SMEX03  (Jackson  et  ah,  2005;  Bosch  et  ah,  2006a). 
The  experiment  was  conducted  over  the  LREW,  Georgia,  and  sites  in  Oklahoma  and 
Alabama.  SMEX03  was  the  second  in  a  chain  of  field  campaigns  conducted  to  validate 
brightness  temperature  data  and  soil  moisture  retrieval  algorithms  for  the  Advanced 
Microwave  Scanning  Radiometer  for  the  Earth  Observing  System  (AMSR-E)  on  the 
Aqua  satellite.  AMSR-E  is  a  passive  microwave  radiometer  sensing  microwave  radia¬ 
tion  on  12  channels  at  six  frequencies  ranging  from  6.9  to  89  GHz  (0.3  to  4.3  cm). 
During  SMEX03,  researchers  collected  aircraft  data  with  the  PSR/CX  sensor  (X-Band, 
2.5-3. 7  cm  wavelength  with  approximately  3  km  spatial  resolution),  satellite  AMSR-E 
X-Band  (0. 3-4.3  cm)  data  (40-75  km  spatial  resolution),  and  in  situ  measured  soil 
moisture  at  19  long-term  sites  and  49  supplemental  sites  (Jackson  et  al.,  2005;  Bosch 
et  al.,  2006a). 

Jackson  et  al.  (2005)  found  that  both  the  PSR/CX  and  AMSR-E  X-band  channels 
were  well  calibrated  and  that  X-band  comparisons  of  the  PSR/CX  high-resolution 
data  with  the  low-resolution  data  from  AMSR-E  indicated  a  linear  scaling  for  the 
range  of  conditions  studied  in  SMEX03.  This  work  was  followed  by  another  analy¬ 
sis  of  AMSR-E  in  the  LREW,  where  LREW  and  three  other  soil  moisture  networks 
were  used  (Jackson  et  al.,  2010).  They  compared  National  Aeronautics  and  Space 
Administration  (NASA)  and  Japan  Aerospace  Exploration  Agency  (JAXA)  soil  mois¬ 
ture  products  to  the  network  observations,  along  with  two  alternative  soil  moisture 
products  developed  using  the  single-channel  algorithm  (SCA)  and  the  land  parameter 
retrieval  model  (LPRM).  They  found  that  the  statistical  results  indicated  that  each 
algorithm  performed  differently  at  each  site.  Furthermore,  the  JAXA  algorithm  per¬ 
formed  better  than  the  NASA  algorithm  under  light-vegetation  conditions,  whereas  the 
NASA  algorithm  was  more  reliable  for  moderate  vegetation.  The  SCA  had  the  lowest 
overall  root  mean  square  error  with  a  small  bias,  while  the  LPRM  had  a  very  large 
overestimation  bias  and  retrieval  errors. 

Sahoo  et  al.  (2008)  estimated  the  soil  moisture  from  AMSR-E  using  the  10.7  GHz 
(2.8  cm)  frequency  channel  and  compared  it  to  the  Land  Surface  Microwave  Emission 
Model  (LSMEM).  In  turn,  they  compared  the  two  methods  to  in  situ  measured  soil 
moisture  datasets  over  the  LREW  and  found  that  the  LSMEM  method  performed  bet¬ 
ter  than  the  current  operational  AMSR-E  retrieval  algorithm. 

DATA 

A  network  of  soil  water  measurement  monitors  has  been  established  in  the 
Suwannee  River  Basin  to  complement  SEWRL’s  30-year  record  of  climatic,  hydro- 
logic,  and  water  quality'  research  (Fig.  2).  Data  are  recorded  from  SEWRL’s  soil  mois¬ 
ture  sensors  every  half  hour  at  three  different  depths:  5.08  cm  (2  inches),  20.32  cm  (8 
inches),  and  30.48  cm  (12  inches)  below  the  surface.  These  data  are  packaged  into  files 
and  processed  for  each  site  (Bosch  et  al.,  2007).  Recorded  information  includes  date 
and  time  of  reading,  soil  temperature,  soil  water  content,  soil  salinity,  and  precipitation 
data.  A  subset  of  these  soil-moisture  network  stations  was  used  for  this  study,  and  the 
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Fig.  2.  Locations  of  soil  moisture  sensors  in  the  LREW. 


number  of  stations  used  was  chosen  to  optimize  hyperspectral  data  collection.  Flights 
across  the  LREW  were  conducted  on  November  30,  2005  and  March  6,  2007.  For  the 
November  30,  2005  flight,  18  soil  moisture  probes  were  chosen,  and  for  the  March  6, 
2007  flight,  29  soil  moisture  probes  were  chosen.8  Flight  lines  were  drawn  over  the 
selected  soil  moisture  meters  for  each  flight  using  United  States  Geological  Survey 
(USGS)  Digital  Ortho  Quads  (DOQ). 

The  aircraft-mounted  SWIR  HyperSpectral  Imaging  (HSI)  sensor,  developed  by 
the  Institute  for  Technology  Development  (ITD)  of  the  Stennis  Space  Center,  was  used 
to  collect  electromagnetic  reflectance  data  over  the  soil  moisture  probe  stations.  The 
system  consists  of  sensors  for  three  different  sections  of  the  electromagnetic  spectrum: 
ultraviolet  (UV),  VNIR,  and  SWIR.  The  sensor  was  mounted  aboard  a  Cessna,  which 
flew  at  a  nominal  6,000  feet  above  the  ground,  resulting  in  a  ground  spatial  resolu¬ 
tion  of  the  image  data  at  just  over  6  meters.  The  extent  of  the  image  in  the  along-track 
direction  was  dictated  by  the  distance  the  aircraft  flew  while  the  sensor  was  recording 
data.  Flight  lines  were  typically  established  by  drawing  multiple  lines  over  a  digital 


"The  first  flight  in  2005  was  used  to  preliminarily  test  the  feasibility  of  the  project.  Based  on  the  results 
from  the  first  flight,  a  more  vigorous  statistical  design  and  analysis  were  conducted  for  the  second  flight  in 
2007. 
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map  of  the  target  area.  The  lines  were  oriented  so  that  all  the  targets  could  be  covered 
within  the  sensor’s  field  of  view  using  the  minimum  number  of  flight  lines.  The  lati¬ 
tude/  longitude  of  the  start  and  stop  points  for  each  line  were  recorded.  The  start/stop 
and  latitude/  longitude  pairings  were  then  uploaded  into  the  GPS  in  the  aircraft. 

The  1TD  SWIR  HSI  sensor  allows  light  to  pass  through  a  lens  and  enter  into  a 
spectrograph  through  a  slit.  The  spectrograph  in  the  sensor  breaks  light  passing  through 
the  slit  into  wavelengths  that  are  in  the  SWIR  range.  The  resulting  light  illuminates  a 
charge-coupled  device  (CCD)  array  and  is  dispersed  with  the  shorter  wavelengths  illu¬ 
minating  the  top  of  the  CCD  array,  and  then  gradually  moving  through  the  spectrum 
with  the  longer  wavelengths  illuminating  the  bottom  of  the  CCD  array. 

The  CCD  array  has  240  lines  and  320  columns.  The  number  of  lines  determines 
the  maximum  number  of  spectral  bands  that  can  be  collected  and  the  number  of  col¬ 
umns  dictates  the  maximum  number  of  cells  of  spatial  information  per  line.  With  this 
SWIR  CCD  array,  240  lines  of  spectral  band  information  could  be  obtained;  however, 
for  this  research  project,  the  number  of  photons  that  are  collected  in  one  I  x  1  cell 
of  the  CCD  array  may  not  result  in  enough  energy  to  excite  the  sensing  element.  To 
accommodate  for  this  possibility,  the  sensor  has  the  capability  to  aggregate  several 
cells  together.  By  increasing  the  number  of  cells  that  are  aggregated  together  dur¬ 
ing  data  collection,  a  larger  area  can  be  created  so  that  the  accumulation  of  photons 
received  by  the  combination  of  cells  can  increase  the  sensitivity  for  each  recording 
element,  thus  improving  the  quality  of  the  data. 

Binning  was  used  to  aggregate  cells  to  create  a  larger  sensing  element  that  records 
the  number  of  detected  photons  for  measurement.  The  two  advantages  to  binning  are: 
(I)  it  increases  the  signal  to  noise  ratio  for  the  data;  and  (2)  binning  the  spatial  dimen¬ 
sion  increases  the  spatial  resolution.  Increasing  the  binning  in  the  spectral  dimension 
decreases  the  number  of  bands  that  is  collected.  For  a  sensor  that  collects  a  maximum 
of  240  spectral  bands,  a  binning  of  the  data  by  two  causes  the  number  of  bands  that 
is  collected  to  be  decreased  by  a  factor  of  two,  leaving  120  bands  of  data.  Likewise, 
binning  the  data  by  four  causes  the  data  to  be  decreased  by  a  factor  of  four,  leaving  60 
bands  of  data. 

For  the  first  flight  on  November  30,  2005,  each  of  the  240  spectral  sensing  lines 
(elements)  was  binned  with  their  neighbors  to  enhance  the  signal  measurement.  This 
resulted  in  every  pair  of  the  spectral  sensing  elements  being  binned  together  for  data 
measurement  and  subsequently  120  bands  of  SWIR  data  were  collected  for  each  flight 
line  in  the  image  data.  By  binning  spectrally  instead  of  spatially,  measurements  of  the 
spectral  absorption  troughs  and  reflection  peaks  were  made  while  maintaining  visual 
identification  of  the  targets  on  the  ground.  High  and  low  bands  outside  the  spectral 
response  of  the  spectrograph  were  omitted  from  the  analysis.  The  resulting  data  set 
contained  85  bands  of  information  with  spectral  response  in  the  range  of  937  to  1700 
nanometers  in  wavelength. 

For  the  second  flight  on  March  6,  2007,  four  adjacent  spectral  sensing  elements 
were  binned  together  to  further  increase  the  signal  measurement.  This  resulted  in 
recording  60,  or  half,  of  the  first  flight’s  number  of  spectral  bands.  All  60  bands  of 
SWIR  data  were  collected  for  each  flight  line  in  the  image  data.  Again,  the  high  and 
low  bands  outside  the  spectral  response  of  the  spectrograph  were  omitted  from  the 
analysis.  Therefore,  the  data  set  that  was  used  for  analysis  contained  43  bands  of 
reflectance  data  with  a  wavelength  range  of  91 0  to  1686  nanometers. 
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Fig.  3.  Calibration  tarps  and  soil  moisture  meter  66  in  the  HS1  image. 


The  data  for  the  second  flight  were  corrected  for  atmospheric  effects  and  system 
noise.  Atmospheric  correction  is  used  to  adjust  for  the  absorption  characteristics  of 
the  atmosphere  and  to  convert  the  data  to  percent  reflectance  (Price,  1987;  Chavez 
and  Mackinnon,  1 994).  It  uses  a  regression  line  with  one  graph  point  representing 
the  raw  digital  number  for  the  dark  point  on  the  x-axis  and  the  corresponding  percent 
reflectance  of  the  rad  iometer  data  on  the  y-axis,  while  the  other  graph  point  represents 
the  raw  digital  number  for  a  bright  point  on  the  x-axis  and  the  corresponding  per¬ 
cent  reflectance  of  the  radiometer  data  on  the  y-axis.  These  regression  equations  were 
used  to  transform  all  the  raw  image  data  to  percent  reflectance.  An  Analytical  Spectral 
Devices  (ASD)  S  W1R  radiometer  was  used  to  capture  radiometric  data  for  comparison 
with  the  hyperspectral  data. 

SWIR  spectra  were  extracted  from  the  radiometer  data  in  percent  reflectance  units 
of  measure  and  paired  with  SWIR  spectra  extracted  from  the  noise-reduced  image 
data  to  generate  a  regression  equation  for  calibration.  The  linear  regression  equation 
from  these  pairings  was  used  to  transform  all  the  image  frames  to  percent  reflectance. 
The  resulting  data  sets  contained  bands  with  a  spectral  response  between  91 0  to  1 700 
nanometers  in  wavelength. 

The  Saviszky  Golay  noise  removal  method  (Press  et  al.,  1997,  p.  994)  was  per¬ 
formed  to  extract  noise  that  is  inherent  in  the  system.  This  method  uses  moments  of 
the  data  values  across  the  spectra  to  maintain  the  signal  while  minimizing  fluctuations 
caused  by  noise  in  the  data,  transforming  the  raw  data  to  a  new  data  set  where  the  noise 
is  minimized.  Calibration  tarps  were  placed  on  the  ground  along  the  flight  line  (Fig. 
3).  The  total  size  of  the  joined  sections  of  tarp  was  36.58  m  (120  feet)  across  track  by 
9. 14  m  (30  feet)  along  track.  Each  section  covered  a  4.57  m  (15  feet)  *  4.57  m  (15 
feet)  area  (thus,  a  total  of  16  sections).  The  tarps  have  eight  reflectance  sections  and 
are  manufactured  to  have  gradually  greater  reflectivity  from  one  section  to  the  next 
(across  track).  They  were  factory  measured  to  have  2%,  4%,  8%,  16%,  32%,  48%, 
64%,  and  88%  reflectance.  However,  in  practice  the  tarp  reflectances  deviate  slightly 
from  these  manufactured  settings.  In  order  to  ensure  a  true  measurement  of  reflectance 
across  each  of  the  tarp  sections,  radiometer  readings  were  taken  on  the  ground  during 
the  time  of  airborne  data  acquisition  (Lewis  and  Finn,  2007). 

For  both  flights,  regions  of  interest  (ROI)  corresponding  with  the  location  of  the 
soil  moisture  probes  were  used  to  extract  the  spectral  information  from  the  SWIR  data 
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Fig.  4.  Average  monthly  precipitation  for  LREW  from  1968  to  2007,  with  2005  and  2007  high¬ 
lighted  (data  from  LREW  database). 


set  for  each  of  the  soil  moisture  probe  locations.  Each  ROI  covers  a  3  x3  pixel  mask 
around  the  soil  moisture  probe  and  the  ground  spatial  resolution  of  the  image  data  is 
just  over  6  meters.  Thus  the  ROI  for  each  soil  moisture  probe  is  approximately  18  by 
18  m.  The  3*3  window  of  pixels  w  ithin  each  ROI  were  then  averaged  to  create  one 
spectral  value  for  each  of  the  soil  moisture  probes.  For  the  soil  moisture  probe  data, 
the  corresponding  time  of  the  HSI  sensor  data  collection  flight  was  used  to  obtain  soil 
moisture  values  for  the  correlative  soil  moisture  probes  used.  These  sensor  data  for  the 
5.08  cm  (2  inches),  20.32  cm  (8  inches),  and  30.48  cm  (12  inches)  probes  were  used 
for  correlation  to  the  SWIR  spectral  information. 

ANALYSIS  AND  RESULTS 

Figure  4  compares  the  average  monthly  precipitation  for  this  area  with  the  monthly 
precipitation  for  the  years  of  data  acquisition,  2005  and  2007.  During  the  months  of 
data  acquisition,  soil  moisture  conditions  were  relatively  dry'  for  November  2005  and 
wet  for  March  2007  (Fig.  5).  Figure  6  shows  an  example  of  the  data  set  collected  over 
the  University  of  Georgia  Tifton  Campus  in  Tifton,  Georgia. 

In  order  to  evaluate  the  utility  of  the  SWIR  as  a  viable  tool  for  soil  moisture 
data  collection,  the  hyperspectral  data  sets  for  the  two  dates  were  used  to  establish  a 
General  Linear  Model  (GLM)  for  each  of  the  three  different  probe  ty  pes.  In  addition, 
an  R 2  analysis  was  used  to  determine  the  strength  of  the  linear  relationship  between  the 
probe  and  SWIR  spectral  band  data.  Recall  from  the  data  section  that  for  the  2005  flight 
each  of  the  240  spectral  sensing  lines  were  binned  with  their  neighbors  to  enhance  the 
signal  measurement.  For  the  2007  flight  four  adjacent  spectral  sensing  elements  were 
binned  together  to  further  increase  the  signal  measurement. 

November  30,  2005 

Data  from  18  of  the  soil  moisture  probes  were  examined  for  the  November  30, 
2005  flight.  Many  of  the  soil  moisture  probes  were  located  in  generally  bare  or  unob¬ 
structed  soil,  but  some  were  nestled  in  significant  ground  cover,  which  made  it  diffi¬ 
cult  to  draw  ROls  with  confidence.  While  none  of  the  soil  moisture  probes  may  have 
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Fig.  6.  SWIR  data  over  the  University  of  Georgia’s  Tifton  Campus  (3I°28'37"  N,  83°3T54" 
W)  located  within  the  watershed,  and  which  served  as  one  of  the  data  collection  areas  (red,  green, 
and  blue  components  are  the  1222,  1106,  and  1563  nanometer  wavelengths,  respectively). 


been  in  100%  bare  soil  and  completely  unobstructed  from  other  ground  artifacts,  six 
soil  moisture  probes  seemed  to  reside  in  significant  ground  cover  or  were  close  to 
other  artifacts  such  as  roadways  or  hay  bales  within  20  meters  of  a  soil  moisture  sta¬ 
tion.  Ultimately,  these  six  entries  of  the  soil  moisture  probe/spectral  data  pairings  were 
removed  from  the  analysis  and  the  remaining  12  pairings  of  data  used  for  the  2005 
analysis. 

Correlations  determined  for  the  12  pairings  of  spectra  per  band  and  the  2-inch 
soil  moisture  data  are  shown  in  Figure  7.  Both  absolute  maxima  and  minima,  and 
local  maxima  and  minima,  can  provide  influencing  locations  across  the  spectra  for 
classification  algorithms.  The  locations  where  the  absolute  and/or  local  maxima  and 
minima  of  the  correlations  provided  points  for  classification  algorithms  w  ere  at  bands 
15,  55,  and  81,  which  correspond  with  wavelengths  of  1062,  1422  and  1664  nanome¬ 
ters  (Table  1 ). 

A  new  GLM  was  generated  to  model  the  5.08  cm  (2  inches),  20.32  cm  (8  inches), 
and  30.48  cm  (12  inches)  soil  moisture  data  with  bands  15  (1062  nm),  55  (1422  nm), 
and  81  (1664  nm).  On  the  2-inch  model,  band  15  did  not  exhibit  a  significance  value 
below  0.05.  Thus,  the  model  was  executed  again  using  only  bands  55  and  81.  The 
results,  shown  in  Table  2,  show  the  2-inch  and  12-inch  models  to  be  significant  at  the 
a  =  0.05  level.  Although  close,  the-8  inch  model  is  not  significant  at  the  0.05  level.  The 
2-inch  model  explained  79%  of  the  variance  in  surface  soil  water  content,  whereas  the 
8-inch  and  12-inch  models  were  less  significant. 
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Correlation  between  spectral  bands  and  2-inch  soil 
moisture  measurement 


Bands 


Fig.  7.  Correlations  between  spectrum  and  2-inch  soil  moisture  from  November  30,  2005 
flight. 


Table  1.  Wavelength  for  Significant  Spectral  Bands 
for  Soil  Moisture  for  November  30,  2005  Flight 


Spectral  band 

Wavelength,  nm 

15 

1062 

55 

1422 

81 

1664 

Table  2.  Model  Coefficients  and  Statistical  Significance  for  Raw  Data 
for  November  30,  2005  Flight 


Model 

R2 

Pr>  F 

inch-2  =  .7728  +  (-.0399  *  b55)  +  (.0107  *  b81) 

0.786 

0.001 

inch-8  =  .3895  +  (.0193  ♦  b55)  +  (.0055  *  b81) 

0.485 

0.050 

inch-12  =  .7549  +  (-.0468  *  b55)  +  (.0148  *  b81) 

0.560 

0.026 

March  6,  2007 

Data  from  29  soil  moisture  sites  and  spectra  pairings  from  the  non-atmospherically 
corrected  data  set  were  used  in  the  GLM  analysis  for  the  2007  data.  The  2-inch  soil 
moisture  data  were  correlated  with  the  spectra  pairings  (Fig.  8). 

Nine  bands  (Bands  29,  30,  32, 33,  34, 36,  38, 39,  and  43  [Table  3])  were  then  used 
to  generate  a  GLM  to  model  for  the  5.08  cm  (2  inches),  20.32  cm  (8  inches),  and  30.48 
cm  (12  inches)  soil  moisture  data.  The  2-inch  and  8-inch  models  are  significant  at  the 
a  =  0.01  level  with  R2  values  above  0.70  (Table  4).  The  12-inch  model  is  significant  at 
the  0.05  level  with  an  R2  of  0.62. 


Correlation 
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Fig.  8.  Correlations  between  spectrum  and  2-inch  soil  moisture  from  March  6,  2007  flight. 


Spectral  bands 


Table  3.  Wavelength  for  Significant  Spectral  Bands 
for  Soil  Moisture  for  March  6,  2007  Flight 


Spectral  band 

Wavelength,  nm 

29 

1299 

30 

1317 

32 

1354 

33 

1372 

34 

1391 

36 

1428 

38 

1465 

39 

1483 

43 

1557 

A  GLM  was  also  generated  for  the  raw  data/atmospheric  calibrated  data  set.  The 
significance  and  R2  are  the  same  (as  in  Table  4),  with  the  model  coefficients  being  dif¬ 
ferent.  These  results  are  shown  in  Table  5. 

CONCLUSIONS 

The  results  of  this  research  show  the  relationship  between  SWIR  data  and  soil 
moisture  and  identify  possible  wavelengths  for  determining  soil  moisture  for  a  given 
soil  type  and  vegetation  distribution.  The  results  for  the  2-inch  model  from  both  flights, 
shown  in  Tables  2  and  4,  were  significant  and  have  a  coefficient  of  determination  ( R 2) 
value  above  0.70.  The  wavelengths  that  best  model  the  soil  moisture  at  the  2-inch 
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Table  4.  Model  Coefficients  and  Statistical  Significance  for  Raw  Data 
of  the  March  6,  2007  Flight 


Model 

R 2 

Pr>  F 

inch-2  =  0.1919  +  (0.007  *  629)  +  (-0.0084  *  630)  +  (0.0015  *  632)  + 
(0.0024  *  633)  +  (-0.0023  *  634)  +  (0.0008  *  636)  +  (-0.0007  *  638)  + 
(00007  *  639)  +  (-0.0008  *  643) 

0.707 

0.0022 

inch-8  =  0.3 102  +  (0.0063  *  629)  +  (-0.0097  *  630)  +  (0.004  *  632)  + 
(0.0027  *  633)  +  (-0.0034  *  634)  +  (0.0008  *  636)  +  (-00006  *  638)  + 
(0.0011  *  639)  +  (-0.0018  *  643) 

0.751 

0.006 

inch- 12  =  0.324  +  (0.0035  *  629)  +  (-0004  *  630)  +  (-0.001 1  *  632)  + 
(0.0035  *  633)  +  (-0.0027  *  634)  +  (0.0013  *  636)  +  (-0.0008  *  638)  + 
(-0.0007  *  639)  +  (-00015  *  643) 

0.623 

0.015 

Table  5.  Model  Coefficients  and  Statistical  Significance  for  Raw  Data/ Atmospheric 
Calibrated  Data  Set 

Model 

R2 

Pr>  F 

inch-2  =  0.1867  +  (2.5693  *  629)  +  (-4.2692  *  630)  +  (1.344  *  632)  + 
(4.001  *  633)  +  (-4.6004  *  634) +  (1.6065  *  636)  +  (-1.4813  *  638)  + 
(1.2715  *  639)  +  (-0.7973  *  643) 

0.707 

0.0022 

inch-8  =  0.3024  +  (2.2851  *  629)  +  (-4.9023  *  630)  +  (3.6488  *  632)  + 
(4.5346  *  633)  +  (-6.8813  *  634)  +  (1.6930  *  636)  +  (-1.1350  *  638)  + 
(2.0093  *  639)  + -1.7591  *  643) 

0.751 

0.006 

inch- 12  =  0.3193  +  (1.2527  *  629)  +  (-2.0254  *  630)  +  (-0.9960  *  632)  + 
(5.7496  *  633)  +  (-5.5 1 58  *  634)  +  (2.8533  *  636)  +  (-1 .6404  *  638)  + 
(1.2887  *  639) +  (-1.3868  *  643) 

0.623 

0.015 

depth  are  in  the  1300  to  1670  range  (Tables  1  and  3).  Of  interest  for  future  work  is  the 
ability  of  particular  bands  to  continuously  predict  soil  moisture  over  time  and  different 
areas,  including  the  strength  and  direction  of  correlation. 

The  data  analysis  illustrates  that  the  ability  for  remotely  sensed  SW1R  data  to  esti¬ 
mate  soil  moisture  deteriorates  at  the  8-inch  and  12-inch  depths.  Although  the  8-inch 
depth  model  had  an  R2  value  of  0.75  and  the  12-inch  depth  model  had  a  slightly  lower 
R2  value  of  0.62  for  the  models  created  from  the  second  flight,  the  first  flight’s  models 
were  not  as  strong  (R2  of  0.49  and  0.56,  respectively).  It  is  not  precisely  known  if 
the  difference  in  R2  values  between  the  two  flights  is  due  solely  to  differences  in  soil 
moisture.  Although  we  suspect  that  is  the  primary  cause,  other  environmental  factors 
differed  between  the  two  flight  dates,  including  atmospheric  conditions  and  vegetation 
cover. 

Ground  cover  can  inhibit  the  ability  for  remotely  sensed  data  to  model  soil  mois¬ 
ture  (e.g.,  Bosch  et  al.,  2006a).  Light  or  the  electromagnetic  wave  needs  to  have 
an  unobstructed  interaction  with  the  soil  and  then  a  clear  path  to  the  sensor  for  the 
SW1R  sensor  to  provide  relevant  data  for  measuring  soil  moisture  at  the  2-inch  depth. 
Logistical  issues  for  future  work  would  include  reducing  the  spatial  resolution  of  the 
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spectral  data  to  minimize  interaction  with  artifacts  such  as  roadways  and  buildings. 
In  addition,  researchers  could  look  more  closely  at  the  relationship  between  reflec¬ 
tance  and  soil  moisture  where  there  is  ground  cover  and  ascertain  the  Normalized 
Difference  Vegetation  Index  (NDVI)  in  those  areas.  This  could  establish  a  minimum 
NDVI  for  ground  cover  that  would  allow  for  estimation  of  soil  moisture.  Alternatively, 
they  could  establish  the  vegetative  (spectral)  response  of  the  ground  cover  to  use  as  a 
proxy  for  soil  moisture. 

The  benefits  of  hyperspectral  data  are  twofold.  First,  the  spatial  resolution  (ca. 
6  m  for  this  study)  is  much  finer  than  the  datasets  used  in  previous  LREW  studies 
by  Jackson  et  al.  (2005)  and  Bosch  et  al.  (2006a)  (ca.  3  km  for  the  PSR/CX  X-band 
and  40-75  km  for  the  AMSR-E  X-band).  In  addition,  it  is  worthwhile  to  note  that 
this  higher  spatial  resolution  (by  an  order  of  magnitude)  is  particularly  important  for 
investigations  of  watersheds,  where  accurate  spatial  representation  can  impact  water 
resource  decisions.  Furthermore,  although  hyperspectral  flights  are  weather  depen¬ 
dent,  as  they  are  mounted  on  aircraft  and  not  satellites,  researchers  have  some  control 
over  the  temporal  resolution.  Therefore,  the  inherent  characteristics  of  the  hyperspec¬ 
tral  data  directly  relate  to  its  strengths  in  determining  soil  moisture. 

Even  though  the  use  of  microwave  sensors  does  warrant  further  research,  the  goal 
of  this  study  was  to  investigate  the  potential  utility  of  hyperspectral  data  with  SWIR 
sensors  to  determine  soil  moisture.  Our  results  indicate  that,  at  least  for  a  2-inch  depth, 
hyperspectral  sensors  can  be  used  to  identify  soil  moisture,  most  notably  in  the  SWIR 
range.  For  situations  where  satellite-mounted  microwave  sensors  are  not  an  optimal 
choice,  especially  in  regard  to  spatial  or  temporal  constraints,  aircraft-mounted  hyper¬ 
spectral  sensors  can  be  a  valid  and  useful  alternative.  Although  the  sensor  used  in  this 
study  was  limited  to  UV,  VN1R,  and  SWIR,  additional  future  studies  using  hyperspec¬ 
tral  sensors  with  microwave  and  thermal  capabilities  are  of  interest,  as  this  study  has 
shown  the  potential  of  hyperspectral  data  for  soil  moisture  evaluation. 

Aside  from  the  microwave  data  sources  mentioned  above  (PSR/CX  and  AMSR-E), 
the  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  can  be  used  to  derive 
the  temperature/vegetation  dryness  index  (TVDI),  which  can  be  used  to  make  infer¬ 
ences  about  soil  moisture.  The  spatial  resolution  of  MODIS  is  finer  than  PSR/CX  and 
AMSR-E  (between  250  and  1000  m),  but  it  is  still  coarser  than  the  hyperspectral  data 
used  in  this  study.  In  addition,  the  temporal  resolution  differences  (satellite  versus  air¬ 
craft)  may  also  be  a  determining  factor  in  sensor  choice.  Although  an  area  of  potential 
future  research  would  be  to  evaluate  the  potential  of  MODIS  to  determine  soil  mois¬ 
ture  in  the  LREW,  it  is  beyond  the  scope  of  this  study. 

The  research  described  herein  provides  evidence  for  the  utility  of  developing 
SWIR  remote  sensing  analysis  to  assess  soil  moisture,  particularly  for  soil  moisture 
measurements  at  the  depth  of  tw  o  inches.  The  ability'  to  use  remote  sensing  techniques 
to  obtain  soil  moisture  information  over  a  large  area  could  decrease  reliance  on  in  situ 
point  measurements,  as  well  as  quickly  provide  data  to  researchers  and  agriculturalists. 
Further  studies  to  test  the  relationships  established  in  this  study  and  to  test  their  applica¬ 
bility  to  other  geographic  areas  could  provide  new  insights  in  soil  moisture  research. 
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